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ABSTRACT 


Lowering the frequency in a superconducting accelerator for a free electron laser (PEL) has the 
potential to reduce the size, cost, and power consumption of the EEL system. A lower frequency 
also enables the use of longer pulses, which has been shown to improve EEE performance. 
Using simulation codes developed at the Naval Postgraduate School, the performance of EEE 
amplifiers and oscillators at several accelerator frequencies is investigated. The results show 
that both EEE amplifier and oscillator performance can be improved by lowering the accelerator 
frequency. In addition, a simulation has been developed that tracks electron trajectories through 
several quadrupole magnets and an EEE undulator. Electron trajectories can directly impact 
EEE performance, and stray electrons may cause damage or harmful radiation if they strike 
beamline components. This simulation has potential as a future research and design tool. 
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CHAPTER 1: 
Introduction 


The free eleetron laser (PEL) is a promising souree of direeted energy (DE) for the United 
States Navy (USN) to eounter emerging and existing threats sueh as antiship missiles (ASMs) 
and unmanned aerial vehieles. A DE weapon delivers energy at the speed of light, eliminating 
the effeetiveness of high-g maneuvers by supersonie missiles. Reeently, the effeetiveness of 
solid state lasers has been demonstrated in a maritime environment at tens of kilowatts, but 
solid state lasers are not expeeted to provide the megawatt power levels that the USN requires 
for ASM defense [1,2]. 

The primary advantages of EEEs over other types of lasers are that they have exeellent optieal 
mode quality, are tunable, and have the potential for good wall plug effieieney [2]. However, 
signifieant ehallenges remain before an EEE weapon system ean be deployed on a ship. One of 
these ehallenges is the overall size and eomplexity of the EEE [1]. 

The eleetron beam used by an EEE is produeed in a linear aeeelerator (linae) powered by radio 
frequeney (RE) radiation. The eurrent in the eleetron beam, the frequeney of the aeeelerator, and 
the eleetron beam quality are all related and impaet EEE performanee. There are engineering 
trade-offs assoeiated with any ehoiee of aeeelerator frequeney; e.g., using a lower accelerator 
frequency has the potential to reduce the size, power requirements, complexity, and cost of 
the linae subsystem [2]. Chapter 3 describes a study I conducted on the effects of accelerator 
frequency on EEE performance. The results of this study will encourage increased consideration 
of lower accelerator frequencies in the design of new EEEs. 

Electron beam generation, acceleration, and transport systems for EEEs are designed to deliver 
high quality electron beams that will result in the best performance. One way to focus the 
electron beam is to use magnetic quadrupoles before the beam enters the section of the EEE 
where the laser light is produced (called the undulator). Once in the undulator, the electrons 
undergo another type of focusing called betatron motion. Chapter 4 describes a simulation that 
I developed to track the trajectories of electrons as they are focused through a set of quadrupole 
magnets and in the undulator. This simulation could be useful for future research or the design 
of focusing systems for a future EEE. 

Eirst, however, I will provide an introduction to EEEs and the theory of their operation. 
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CHAPTER 2: 

Free Electron Laser Theory 


2.1 What is an FEL? 

A free electron laser (FEL) is a laser which has a gain medium composed of free electrons— 
electrons that are not bound in atoms or molecules. This distinguishes FELs from the thousands 
of other laser designs and gives FELs some unique properties and advantages not found in 
conventional lasers. 

2.1.1 Conventional Lasers 

The term laser is an acronym for Light Amplification by Stimulated Emission of Radiation. 
Stimulated emission is a quantum-mechanical process in which excited electrons emit photons 
with a certain energy and phase in the presence of other photons with the same energy and 
phase. Stimulated emission is closely related to another quantum-mechanical process called 
spontaneous emission, which occurs when an electron transitions to a lower energy level spon¬ 
taneously and emits a photon with a random phase [3]. 

A simplified conventional laser consists of an energy source, a gain medium, and an optical 
resonator. The energy source excites the electrons in the gain medium to begin the process of 
spontaneous emission. The resulting photons resonate in the optical cavity, providing feedback 
and driving stimulated emission, amplifying coherent electromagnetic radiation. Usually, one 
side of the optical resonator will be partially transmissive to outcouple the laser beam [3]. 

The gain medium determines the optical properties of the laser, such as its wavelength, and it 
also might limit the amount of power that can be produced before the gain medium is damaged. 

2.1.2 FEL Attributes 

An FEL, on the other hand, can be designed for any wavelength from microwaves to x-rays and 
is not power-limited by damage to its gain medium, which is an electron beam that is continually 
replenished at the speed of light. FELs are also more flexible than conventional lasers because a 
single design is tunable over a larger range of wavelengths than any other laser; a conventional 
laser has a fixed wavelength or small range of wavelengths determined by the energy levels of 
the gain medium [1]. 
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The primary disadvantages of FELs are their size and eost. Laboratory FELs are typieally 
several meters long, depending on their wavelength and eonfiguration. Cost will vary depending 
on applieation, desired output power, and other faetors, but the eost per engagement is low (the 
eost of a few gallons of ship fuel, for example) [4]. 

2.1.3 Naval Applications 

The primary naval applieation of an FEE is antiship missile defense (ASMD). Antiship missiles 
(ASMs) pose a signifieant threat to the surfaee navy, and as they beeome faster, harder to deteet, 
and inereasingly deadly, a similarly advaneed and eapable defensive weapon is required [1]. 

FELs deliver lethal energy at the speed of light—almost instantaneously—eliminating the ad¬ 
vantage gained by maneuvering to evade defensive measures. Many ASMs evade eurrent de¬ 
fensive systems by performing high “g” maneuvers in their terminal phase; sueh maneuvers are 
easily eountered by a system that delivers its ordinanee at the speed of light [1]. 

FELs are able to destroy ineoming missiles or aireraft quiekly and surgieally. Supersonie ASMs 
ean eover the distanee from the horizon to their target in less than a minute. A megawatt-elass 
FEE will aehieve an aetive kill in a matter of seeonds, allowing the weapon to quiekly engage 
other threats. Targeting and traeking speeds are limited only by the radar and direetor systems 
in use, allowing rapid target switehing and surgieal delivery of lethal energy. 

An FEE might also be used against a target requiring a graduated response, sueh as an ineoming 
small eraft. Low power levels eould be used to deliver a warning, with inereased energy levels 
disabling the eraft or detonating ordinanee on a suieide eraft, for example [4]. 

2.1.4 Major Components 

An FEE eonsists of three primary eomponents (shown in Figure 2.1): an injeetor, where the 
eleetron beam is produeed at a eathode; a linear aeeelerator (linae) where the eleetron beam is 
aeeelerated to highly relativistie energies; and an undulator, where the FEE interaetion produees 
laser light [2]. 


Injector. The injeetor eonsists of the eathode, whieh produees the free eleetrons, and an aeeel- 
erating eavity whieh quiekly aeeelerates the eleetrons up to relativistie speeds (~5MeV) and 
into the linae [2]. 
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Figure 2.1: An FEL consists of three major eomponents: an injeetor, a linae, and an undulator. 
An FEL oseillator has an optieal eavity ereated by two mirrors, one of whieh outeouples a small 
amount of the light. In an FEL amplifier, a seed laser is used to generate the initial laser beam 
and there is no optieal eavity. After [4]. 


Linac. The linae further aeeelerates the eleetrons until they are highly relativistie (~100MeV 
for an infrared weapons-elass FEL). The linae usually eonsists of a series of radio frequeney 
(RE) eavities through whieh the eleetrons are eontinuously aeeelerated by the alternating eleetrie 
fields. In order to inerease the effieieney and gradient of the RE eavities, supereondueting 
materials are often used to eonstruet the eavity walls, and a eryoplant is required to eool the 
eavity walls to supereondueting temperatures (2 K to 4 K in this ease) [2]. 


Undulator. An undulator (also ealled a “wiggler”) is made up of permanent magnets or elee- 
tromagnets that ereate a transverse periodie magnetie field along its length. The transverse 
magnetie fields eause eleetrons passing through it to oseillate and emit radiation. The light that 
is emitted in this proeess has a mueh shorter wavelength than the period of the magnetie field in 
the undulator [2]. 


2.2 Types of FELs 

There are two primary elasses of FELs: oseillators and amplifiers. Eaeh has its own advantages 
and disadvantages, and both FEL oseillators and FEL amplifiers are in use at researeh faeilities 
throughout the world [5]. 
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2.2.1 Oscillators 

An FEL oscillator has an undulator within an optical cavity. The light in the optical cavity is 
repeatedly amplified by the FEL gain process over many passes. In order to outcouple a portion 
of the light, one of the cavity mirrors is typically partially transparent [2]. 

FEL oscillators have a few notable advantages over amplifiers. Perhaps most important for 
naval application is that they require significantly shorter undulators than amplifiers—often by 
an order of magnitude for the same output power [2]. Oscillators have superior optical mode 
quality at their output, which improves beam propagation through the atmosphere to a target 
[1]. Simulations have shown that oscillators tolerate vibrations better, which is a significant 
advantage for a shipboard application [6]. Finally, oscillators are less sensitive to electron beam 
quality due to their shorter undulators. 

There are a few notable disadvantages of oscillators related to their optical cavities. The first 
is that a suitable material must be found for the optical cavity mirrors that can tolerate high 
power radiation at the wavelength of the FEL; this is especially true for the outcoupling mirror. 
Second, it is typically necessary to provide cooling to the mirrors at high optical power [2]. The 
third disadvantage is that the choice of mirror material might limit the tunability of the FEL if 
the reflectivity of the material varies significantly in the range of wavelengths of interest [1]. 

2.2.2 Amplifiers 

FEE amplifiers have no optical cavity; instead, a seed laser is used to inject an initial optical 
beam. The seed laser fixes the wavelength of the FEL, making it untunable (unless the seed laser 
is changed). A longer undulator is also required in order to achieve significant amplification of 
the optical field over a single pass [2]. 

Without an optical cavity, an amplifier has the advantage that there is no outcoupling mirror to 
damage at high power. However, an optical element must be used at some point nearby to direct 
the output beam, so amplifiers actually share many of the disadvantages related to the optical 
cavity of an oscillator [2]. 

2.3 FEL Theory 

2.3.1 Electron Equations of Motion 

We will begin our foray into the classical FEL theory by deriving the equations of motion 
and the energy equation for the electrons in the undulator. It is assumed that the reader has a 
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basic understanding of relativity and classical electromagnetism. The relativistic Lorentz force 
equation for a particle with charge q and mass m is^ 

= —(e + JSxb), ( 2 . 1 ) 

at me V / 

where 7 is the Lorentz factor, P = v/c is the particle’s relativistic velocity vector, c is the speed 
of light, E is the electric field vector, and B is the magnetic field vector. The work done on an 
electron by the fields of the undulator and laser light is given by 

r=—E-i3. ( 2 . 2 ) 

me 

To take these equations any further, we will need to specify the electric and magnetic fields 
inside the undulator, which are a result of both the undulator fields and the light. The most 
commonly used undulator in practice is a linear undulator in which a series of alternating pole 
faces create a sinusoidal and transverse variation in the magnetic field along the undulator. 
However, the field created by a helical undulator is much easier to deal with analytically. A 
helical undulator creates a magnetic field that spirals along the undulator and has the form^ 

B„ = Bq [cos(koz)x + sin(koz)y] (2.3) 

where Ao = 2 ;r/ko is the period of the undulator field (usually a few cm), and z is the longitudinal 
position in the undulator [2]. The light produced in such a helical undulator will be helically 
polarized and will have electric and magnetic fields 

E = E (cos V/x — sin t//y), (2.4a) 

B; = E (sint//x + cos V/y), (2.4b) 

where = kz — cot and k = 2n/X. The terms X and to are the usual optical wavelength and 
angular frequency, respectively. 

We can now insert Equations (2.3), (2.4a), and (2.4b) into Equation (2.1). After a bit of algebra 
(applying the assumption that ^ 1 ), we arrive at equations for the transverse and longitudinal 

'This is the Lorentz force equation in Gaussian units (which are used throughout this thesis). For the reader un¬ 
familiar with this system of units, the classic graduate electromagnetism text by Jackson [7] provides a reasonable 
primer. 

^We will orient the undulator along the z axis with the electron bunches traveling in the positive z direction. 
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motion of the electrons: 


dirP±) 

dt 


djrl^z) 

dt 


e 

me L 
e 

me 


-l^zBosm{koz)x + PzBocos{koz)y , 

I Px [5o sin(fcoz) + E cos i//] - Py [5o cos(fcoz) + £ sin t//] |, 


(2.5a) 

(2.5b) 


where we defined = (j8;c,j8y) and j3 = We also substituted q = —e, where e is 

the elementary electron charge magnitude. Integrate Equation (2.5a) by inspection to get the 
transverse electron velocity vector, 


P 


K 

y 


cos(fcoz)x + sin(fcoz)y], 


( 2 . 6 ) 


where K = eBoX^/lnme^ is the dimensionless undulator parameter. For most FELs, K is near 
unity [2]. By inserting Equations (2.4a) and (2.6) into Equation (2.2) we get the rate of change 
of the electron energy in the undulator of an FEE: 

eKE 

7 ^-cos(^ + ^), (2.7) 

yme 

where we defined the electron phase, ^ = {kQ-\- k)z — cot, and (j) is the optical phase. The electron 
phase describes the position of the electron on the scale of an optical wavelength. Note that 


= {ko + k)z - cot 


/In 2n\ 

\ Ao A / 


z—cot 



z —cot 


because Aq 3> A. The term ^ thus has two significant terms: {2n/X)z, which gives the phase of 
the electron relative to an optical wavelength, and cot, which accounts for the changing electric 
field of the optical wave. 


2.3.2 The FEL Resonance Condition 

We are now poised to derive one of the most important equations in FEE physics: the resonance 
condition. The FEE resonance condition gives the approximate wavelength of the light an FEE 
will produce. 

Electron-Photon Race 

Figure 2.2 shows the path of an electron as it oscillates through one period of the undulator. The 
electron will emit roughly one wavelength of light over each of these periods, and because the 
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Figure 2.2: The electron-photon race. As the electron passes through roughly one undulator 
period it emits one wavelength of light. This combined with the difference between the electron 
and photon velocity leads to the resonance condition. 


electron is traveling slightly slower than the light, it slips back by a distance A relative to the 
light over each undulator period. 

The difference between the photon and electron velocity is c(l —/3^) and the duration of the 
race over one undulator period is At = Xq/^^c. The winning distance is then 


A 


c(l-/3z)At 


or 

A=Ao^^^, (2.8) 

Pz 

which is the PEL resonance condition. This expression is exact, but a simple approximation 
results in a more useful form. Recall that 7=1/ and = [5^ + Py + P^ so we can 

write 

2 = 1 - = 1 - (A.- + /32 + ft^) = = 

where we used Equation (2.6) to get the squared magnitude of Solving for P^, we get 
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(applying a binomial expansion with the assumption that 73> 1 ) 


#5= = 



yl 


(l+i^2) 


(2.9) 


Substituting Equation (2.9) into Equation (2.8) we have 


A 


Ao 

Ao 


(1-/3J 



I'f- 


1 


or 


A Ao 


(l+i^2) 

If 


( 2 . 10 ) 


Erom Equation (2.10) we see that the eleetron beam energy {ymf) and the undulator parame¬ 
ters (Ao, K) determine the wavelength of an EEE. The undulator parameter, K, is usually near 
unity and sinee 7 ^ 1, it follows that the wavelength of the light produeed by an EEE will be 
signifieantly shorter than the undulator period, Ao [2]. The resonanee eondition also highlights 
the tunability of EEEs: it is relatively easy to ehange 7 (and in some eases K) to produee a range 
of wavelengths with a single deviee, something traditional lasers eannot do [ 8 ]. 


2.3.3 The FEL Pendulum Equation 

The EEE pendulum equation deseribes how eleetrons evolve in phase spaee. In order to simplify 
our work, we define a few dimensionless quantities. 

We define the dimensionless “time,” T, so that x = ct/L where L is the length of the undulator. 
As an eleetron traverses the undulator, x goes from t = 0 to T = 1. We make use of a speeial 
notation for a derivative with respeet to dimensionless time; an open eirele is used instead of a 

O 

solid dot, like (■ ■ ■) = d{ - ■ ■)/dx. 

The eleetron phase veloeity, v, is the derivative of the eleetron phase with respeet to dimension¬ 
less time and is 

^ ^ ^ ^ ^ ^ =^[(^0 + ^)^-^], (2.11) 
sinee z = /3^c and k = (o/c. We ean take another derivative with respeet to dimensionless time 
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and get 


= + (2.12) 

c at c \ / c 

because k 3> ko. Taking a time derivative of Equation (2.9), we get 

(2.13) 

Near resonance, we can write 

r 2/2 

h 9 TT _ 

By plugging Equations (2.13) and (2.14) into Equation (2.12) we get 

o L r 

V ^ -4;rA^-. 
c 7 

Substituting Equation (2.7) into Equation (2.15), gives the EEE pendulum equation [9] 

V ^ |a|cos(^+ ^), (2.16) 


(2.14) 


(2.15) 


where we defined the dimensionless optical field 

, , 4;rA^ei^L|E| 

1^1 2 7 ’ 
y^mc^ 


(2.17) 


where we replaced y with yo, assuming that the electrons do not lose much energy over the 
length of the undulator. 


Exploring the Pendulum Equation 

The pendulum equation describes how electrons in the undulator evolve over time in dimen¬ 
sionless phase space, and helps in gaining an intuitive sense of the EEE interaction between the 
light and the electrons. 

Electron phase. The electron phase, is a measure of the position of the electron with respect 
to the optical wave on the scale of an optical wavelength, A. A real electron bunch has on the 
order of 10^ electrons that are initially randomly positioned in each optical-wavelength-thick 
“slice.” 
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(a) Small closed orbit. (b) Large closed orbit. (c) Open orbit. 

Figure 2.3: Phase space plots for the motion of a pendulum. For small amplitude oscillations, 
the phase space orbit will be closed as in (a). If the pendulum swings with amplitude great 
enough to reach the top but not go over, the phase space orbit will still be closed, and will look 
like (b). If the pendulum swings up and over to the top, its phase space orbit will be open, as in 
(c). 


Phase velocity. The rate of change of the electron phase, v, is 

v = l=L[iko + k)P,-k], 


and from Equation (2.15) we can approximate 

Ay 

Av 4nN —. 
7 


(2.18) 


Note that the change in electron phase velocity is proportional to a change in electron energy. 
To transfer energy from the electrons to the optical field, the average energy of the electrons^ 
must decrease, which corresponds to a requirement that (Av) < 0. 


Phase space trajectories. There are two types of trajectories in phase space: open and closed. 
Consider a classical pendulum with the well known equation of motion 

0 ~ —y sin0, 

where g is the acceleration due to gravity and I is the length of the pendulum. For relatively 
small oscillation angles, the pendulum will have a closed trajectory in phase space. If the 
pendulum is swinging with enough energy to continue swinging over the top of its anchor, 
however, it will have an open orbit. Examples of each of these cases are shown in Figure 2.3 [4]. 

^We will use the notation (• • •) to indicate an average over a number of sample electrons throughout this thesis. 
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V 



Figure 2.4: The separatrix for the eleetron phase spaee of an PEL separates the open and elosed 
phase spaee orbits. The separatrix is shown as a blue dashed line, and the eleetron orbits are 
shown in red. 


It is useful to define the separatrix, which separates the open and closed orbits in phase space 
for a particular system. Initial system conditions determine whether the pendulum (or electron) 
is inside or outside the separatrix with a closed or open phase space path, respectively. An FEE 
has a separatrix [8] defined by 


V^^ = 2|a| [l+sin(^^+ <)>)], (2.19) 

where the subscript “s” indicates that these values of v and ^ lie along the separatrix. An FEE 
separatrix with open and closed phase space orbits is plotted in Figure 2.4. 

We note that the separatrix depends on the dimensionless optical field, \a\, and the optical phase, 
(j). These values will change in a typical FEE over time and as the electrons travel through the 
undulator, so the separatrix is not fixed. In fact, it can be changed intentionally to improve the 
performance of some FEEs, a technique we will discuss in Section 2.4.3. 

Phase space trajectories of multiple electrons. Now that we know how individual electrons 
evolve in phase space, we consider how the ~10^ electrons in each A-thick optical slice behave. 
An FEE needs two things to lase: the electrons must lose energy to the optical field on average 
(giving power to the optical field) and they must bunch on the scale of the optical wavelength 
(emitting coherent radiation) [2]. 
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(a) Electron evolution at resonance. (b) Electron evolution above resonance. 

Figure 2.5: Electron evolution in phase space. When electrons enter the undulator at resonance 
(Vq = 0), the EEL will have no gain because half of the electrons will lose energy while the other 
half gain energy as in (a). If electrons are injected into the undulator slightly above resonance 
as in (b), where Vq = y/a, they bunch and lose energy on average (Av/v —0.04 in this case). 
To produce these phase space plots, 21 sample electrons were allowed to evolve according to 
the pendulum equation, Equation (2.16). Initial (t = 0) and final (t = 1) electron positions are 
shown in green and blue, respectively. The electron paths are shown in a gradient from yellow 
to red, and the separatrix is shown as a dashed black line. 


Consider an EEE in which all of the electrons enter the undulator on resonance (Vq = 0). As 
shown in Eigure 2.5a, though the electrons bunch on the scale of an optical wavelength and 
have maximum energy exchange, there is virtually no net energy exchange and the EEE will 
not work. 

However, if the electrons are injected slightly above resonance (Vq > 0), as shown in Eig¬ 
ure 2.5b, the electrons still bunch as desired but now lose energy on average. Injecting the 
electrons below resonance (Vq < 0) would also cause the electrons to bunch, but they would 
gain energy from the optical field instead of giving their energy to it. 


2.3.4 Optical Field Evolution 

We have now described in some detail the motion of electrons in the undulator of an EEE and 
how they exchange energy with the optical field. We are missing the equations governing the 
evolution of the optical field in an EEE—the EEE wave equation. With it, we will have a more 
complete picture of the EEE interaction. 
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Low Gain in Weak Optical Fields 

Before deriving the FEL wave equation, we eonsider a speeial ease: low gain in weak optieal 
fields. From eonservation of energy, we know that the energy lost by the eleetrons in the undu- 
lator is gained by the optieal field. To find the ehange in the eleetron energy, we iteratively solve 


C = J vdt 

°° o 

^ = V ^ |a|eos(^ + ^)- 


A seeond-order solution in \a\ is required, and the result is derived in [9] as 


(Av) 


|«o|" 


eos(VoT) - 1 + -VoTsin(voT) 


( 2 . 20 ) 


whieh is valid for low gain (when a ^ao does not ehange mueh) and weak fields (ao 7t). The 
total eleetron energy ehange per buneh is 


average energy lost per eleetron 

X 

number of eleetrons 

= 

(Ay) mc^ 

X 

pV 


where p is the eleetron density and V is the volume oeeupied by the buneh. The energy in the 
optieal field is E^V jAn and the gain is the total energy ehange divided by the energy in the 
optieal field: 


G = 


-r(Av) 2 
-me 


AkN 


(pv) 


Ak 

Wv 


—ymcp 

NE^ 


(Av), 


where we used Equation (2.18) to approximate (Ay) and inserted a minus sign to make the gain 
positive when the optieal field grows (when Av is negative). Inserting Equation (2.20) into the 
above expression gives 


^ 2 Ne 2 “ 2eos(VoT) - VoTsin(voT)] 


''0 


yme^p \6n^N^e^K^L^E^ 
2NE^ Vofi 


me 


2\2 


2 - 2eos(VoT) - VoTsin(voT)] 


SK^Ne^K^L^p 
VQy^me^ 


2 - 2eos(VoT) - VoTsin(voT)] 
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Figure 2.6: Gain curve in weak fields (\a\ < k) with low gain {j < k). Optieal gain (G) is 
plotted versus average initial eleetron phase veloeity (Vq). The gain eurve is antisymmetrie 
about Vo = 0 and has a maximum around vq ~ 2 . 6 . 



where we replaeed 7 with 70 , the initial eleetron energy, sinee we are assuming low gain (the 
eleetrons lose very little energy over the length of the undulator). Define the dimensionless 
current density 


STt^Ne^K^L^p 

Y^mc^ 


( 2 . 21 ) 


and we have 


G 



2eos(voT) - VoTsin(voT) 


( 2 . 22 ) 


The plot of Equation (2.22) in Figure 2.6 reveals what we knew already: there should be no 
gain if all of the eleetrons are at resonanee (Vq = 0). Note that the gain curve is antisymmetrie 
and that the maximum gain oeeurs at Vq ~ 2.6 [4]. 


Development of Coherence 

Coherent light in an FEE is produeed by eleetrons radiating in phase and at the same wavelength. 
The development of eoherenee in FEEs requires that the eleetrons buneh on the seale of an 
optieal wavelength and lose energy (on average) in order to transfer energy to the optieal field 
[ 2 ]. 

In an FEE oseillator, the initial field is broadband radiation from spontaneous emission, and 
eoherenee develops by mode competition over many passes through the undulator [10]. For 
example, in weak fields and with low gain, the strong peak in the gain eurve plotted in Figure 2.6 
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means electrons at that energy would lose more average energy than those at lower or higher 
energy,^ ensuring that over many passes, a strong peak would develop in the corresponding 
optical spectrum. 

Unlike an oscillator, an amplifier does not have an optical cavity that allows coherence to de¬ 
velop over many passes [9]. An amplifier is seeded with an optical pulse that has a finite spread 
in wavelengths, and this pulse is critical in determining the output wavelength and coherence. 
However, amplifiers are not reliant on spontaneous emission, and with quality seed pulses, they 
are capable of producing highly coherent pulses [11]. 

Ultimately, the coherence of a pulse will be limited by its length, which is related to the electron 
pulse length and the slippage distance, NX. Fourier analysis gives a fundamental limit on 
the frequency content of a pulse with a given length which requires that Atrms^OJrms ~ 1/2, 
where Atrms is the root mean square (rms) pulse duration and AcOrms is the rms spectral width. 
When the spectrum of a pulse is as narrow as possible, the pulse is said to be “transform- 
limited” [1,11,12]. 

FEL Wave Equation 

In Section 2.3.4 we developed an expression for low gain in weak fields. For more general cases, 
a wave equation is required to describe the evolution of the optical field. As is common in laser 
physics, we make the assumption that the electric field envelope is slowly varying in space on 
the scale of its wavelength and slowly varying in time compared to the optical frequency. This is 
called the “slowly-varying envelope” approximation [9]. We begin with Maxwell’s full vector 
wave equation in the Coulomb gauge,^ 

= ( 2 . 23 ) 

where V is the Laplacian operator and A(x, t) is the vector potential from which we can recover 
the electric and magnetic fields using 

E = -i^ 

c dt 

B = VxA. 

“^The electron energy is closely related to its phase velocity, v. 

^Again, in Gaussian units. See [7] for information on this system of units. 
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The term on the right side of Equation (2.23) is the transverse eurrent density, 


n 


where is the veloeity of the eleetron at position r„ and 5*^^^ (■ ■ ■) is the Dirae delta funetion. 
We will expand the eurrent density term (J^) using the equations of motion derived for the 
eleetrons in Seetion 2.3.1. The eurrent density is 

n 

= -ec£ 

n 

= 

n 

= ec^Re|/e^——r„) 


_^ 

- {cos(koz)x + sin(fcoz)y) 

, 7 


5(3)(x-r„) 


ke 


-K 

y 




\5W( 




x-r„ 


or 


^ ecKpK& 



(2.24) 


where £ = —ix + y is the eomplex polarization veetor of a eireularly polarized souree (sueh as 
that produeed by a helieal undulator), we used Equation (2.6) and we defined p as the number 
of eleetrons per unit volume for a small volume element dV at (x,t). The volume element 
dV is large enough to eontain a signifieant number of eleetrons but small enough to span only 
a few optieal wavelengths where the field envelope is eonstant. The summation is over the 
n eleetrons in the volume element. In Equation (2.24) we made an approximation where we 
averaged over all of the eleetrons in the volume element dV ; this is eonsistent with the slowly 
varying amplitude and phase approximation we are using for the laser light. 


Eor now, we ignore diffraetion (we diseuss diffraetion in Seetion 2.4.5) and assume the follow¬ 
ing eomplex solution to Equation (2.23): 

A(x,t) ^ A(z,t) ^ 

where we made the approximation that the solution is independent of the transverse eoordinates 
sinee we are ignoring diffraetion. The Eaplaeian term on the left hand side of Equation (2.23) 
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simplifies considerably now that A = A(z, ?) is not a function of x or y. We define 



i{kz—cot) 


SO that A = and we have 


dz 


d^A 

dz^ 


-z—^*£+E^—e 
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^+ikE^ 
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V dz 
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where we dropped the second-order spatial derivative term because it is small compared to the 
first derivative term with the slowly varying envelope assumption. Similarly, the temporal term 
is 


dA 

d^A 

~W 

d^A 

~W 


fdE 

-^-icoE 1 'Pe 
V dt 


dE 


2 d^E 

-2i(0^ - (o E + ^^ 

dt dt^ 

-2ico^ - cofE^ 'Pe, 




where we have dropped the second-order time derivative term since it is small compared to the 
first derivative term with the slowly varying envelope assumption. Inserting the spatial and time 
derivatives of A(z,t) as well as the current density from Equation (2.24) into Equation (2.23), 
we get 


dE 2 
2ik^ - kE -f 


2ico dE (ef 


dz dt 


Ak 


H— ttE =- ecKpie ( 


-ikcfZ 


2i (^ -AneKpie ( 


dz c dt 


\ r 

—ikoz ’ 


\ r 


where we used the complex forms of A(z,t) and J^, recognizing that physical quantities are 
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represented by real parts of each. The polarization vector £ appears on both sides, so all of the 
components on each side must be equal. If we multiply both sides by and drop the 

polarization vectors, we get 


2k 




E — —AnekKp 


y-i[{k+kQ)z-(Ot] 


7 


(2.25) 


where we brought the term inside the average because we agreed that the volume 

element dV is small enough to only span a few optical wavelengths, and the optical parameters 
are slowly varying on this scale. We simplify Equation (2.25) further using a technique called 
the “method of characteristics” where we introduce a coordinate u = z — ct which follows the 
light [9]. This allows us to eliminate the partial derivative with respect to z and Equation (2.25) 
simplifies to 


2k dE 
c dt 


= -AnekKp 


-i[{k+kQ)z-(M] 


7 


(2.26) 


which describes how the electric field changes with time in response to the distribution of 
charge. We previously defined ^ = {kA-kQ)z — cot and the dimensionless time T = ct/L. We 
substitute and get 


2kcdE , ^ 

-=— = -AnekKp ( — 

cL dx 


K' 


\ ^ 


dE _ ! e 


-K' 


^ = -2neKLp 
dx \ 7 


Multiplying both sides by AnNeKL/y^mc^. 


d (AnNeKLE\ 
dx\ YqIuc^ ) 

a 


(%Tl^Ne^K^E^p\ 
V y^mc^ J 

'' -V-' 

fro 



So 


O 



(2.27) 


where (■■■) = (9(---)/(9 t is a partial derivative with respect to dimensionless time, 7o is the 
initial electron energy, a is the dimensionless optical field and j is the dimensionless current 
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density. If we make the assumption that the electrons do not lose much of their energy as they 
travel through the undulator ( 7 o ~ 7 ) we can simplify Equation (2.27) to 



which is the EEL wave equation [9]. 


(2.28) 


2.3.5 Understanding the FEL Wave Equation 

The EEE wave equation reveals the interaction between the optical field and the electrons in 
the undulator. It shows how the dimensionless optical field changes based on the dimensionless 
current density and the bunching term, We will now examine the effects of each of 

these terms in detail. 


Dimensionless Current Density 

The dimensionless current density was defined in Equation (2.21) as 

J — Q 9 • 

Y^mc^ 

Erom Equation (2.28) we see that the larger j is, the faster a changes. The dimensionless current 
density has the following dependencies: 

• i°^K^ 

. jocL^oc NL^ (since L = NAq) 

• jocp 

These dependencies have significant ramifications in the design of EEEs. Recall that from the 
resonance condition. Equation (2.10), A 7 ^^. Increasing the electron beam energy (ymc^) is 
the primary means for shortening the wavelength of an EEE [13]. Doing so, however, comes 
at a significant cost to j, since the dimensionless current density is proportional to 7 ^^. The 
dimensionless current density must be large for a strong interaction to occur, so in most cases 
it is necessary to increase L, the length of the undulator, to compensate for increases in 7 . Eor 
example, the SPRING -8 x-ray EEE in Japan has a wavelength of ~1 A and an undulator section 
that is ~ 100 m in length [5]. 
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Dimensionless Field 

The dimensionless optieal field was defined in Equation (2.17) as 

AnNeKLE 

U = 2 2 ■ 

7Qmc^ 

Like the eomplex eleetrie field E = \E\e''^, ah also eomplex, and has an amplitude and phase: 

a = Ur + iai = \a\e“^j 

where Ur and ai are the real and imaginary parts of a, respeetively; (j) = tan ^ (a,-/a^) is the phase 
of a (and E); and \a\ = \JalA-a^. The rate of ehange of a is 

a = \a\e'^ + \a\e'^i(^ ~ — j i 

where the last equality is from Equation (2.28). Rearranging, we have 

|a| + \a\ — 

= -j (cos(C + (^))) +/j (sin(C + 0)). 


By equating the real and imaginary parts we get 


a\ 

O 

(l> 


-j{cos{C+ <!))) 


(2.29a) 

(2.29b) 


Bunching Term 

Assuming that the optieal 
affeets the optieal 
positioned, and 


phase ^ 0 for simplieity, we investigate how the bunehing term 

gain. At the start of the undulator, the eleetrons will be randomly 



so there will be essentially no ehange in a and no gain. 


As the eleetrons pass through the undulator they will begin to buneh aeeording to the pendulum 
equation. Equation (2.16), and the bunehing term will be nonzero. If the eleetrons buneh near 
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C = ^, then from Equation (2.29a) we see that 


\a\ — —jcosn = j, 

which implies that the optical field will grow [2]. On the other hand, if the electrons bunch near 

O 

C = 0, then \a\ ~ —j and the optical field amplitude will decrease. Similarly, if the electrons 
bunch near ^ = rLKjl there will be no change in the amplitude of the optical field; instead, the 
optical phase will evolve according to Equation (2.29b) [14]. This result shows what was stated 
in Section 2.3.3 without proof: electrons must bunch on the scale of an optical wavelength and 
lose energy for an EEE to work. 

FEE Wave Equation and the Pendulum Equation 

The pendulum equation. Equation (2.16), and the EEE wave equation. Equation (2.28), form a 
pair of coupled differential equations that describe the feedback that occurs between the elec¬ 
trons and the light in the undulator. The pendulum equation describes how the amplitude of the 
optical field drives the evolution of the electrons in phase space, and the EEE wave equation 
reveals how the optical phase and amplitude will evolve in response to the distribution of the 
electrons in phase space. 

2.4 Practical Design Considerations 

When designing an EEE, issues can arise that are not addressed by the basic theory we have 
discussed up to this point. In this section we consider some of these practical design concepts, 
and in a few cases, discuss methods to address them. 

2.4.1 Electron Beam Quality 

The development of high quality, high current (~1 A) electron beams is an area of active re¬ 
search for EEEs. Major challenges exist in the design of the cathode, injector, and accelerator 
systems involving component lifetime, cooling, structural integrity and vibration tolerance [2]. 

The quality of the electron beam can affect EEE performance considerably. The evolution of the 
optical field depends strongly on the electron density, as we showed in Section 2.3.5; as such, 
poorly formed electron bunches result in degraded performance. Similarly, spectral quality of 
the optical pulse depends on the energy spread of the electrons. Electron beam quality will be 
discussed more carefully in Chapter 3. 
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2.4.2 Betatron Motion 

Electrons passing through an undulator undergo a slower betatron motion in addition to the fast 
“wiggling” motion caused by the transverse periodic magnetic field in the undulator. This beta¬ 
tron motion is caused by the transverse variation of the magnetic field strength in the undulator. 
The equation of motion is that of a simple harmonic oscillator and has the form [15] 

oo fKkolV 2 

where we define the dimensionless betatron frequency as 

KkoL IkNK 

COa = - = -. 

P Y Y 

Typical betatron frequency values are on the order of 111, corresponding to one oscillation over 
the length of the undulator (t = 0 —)■ 1). We will discuss betatron motion in more detail in 
Chapter 4. 


2.4.3 Tapered Undulators 

In strong fields (|a| > ;r), it is possible for the gain to saturate due to over-bunching of electrons. 
When \a\ is large, many electrons can become trapped in the separatrix^ and begin to take 
back energy from the optical field. Left unchecked, the electrons will evolve according to the 
pendulum equation, executing synchrotron oscillations and alternating between giving energy 
to the light and taking it away [8,9]. 

To prevent over-bunching, we can alter the resonance condition along the undulator, which 
results in a new form of the pendulum equation and changes the shape of the separatrix. The 
most straightforward method of achieving this is to change the undulator parameter {K) along 
the length of the undulator by decreasing the magnetic field on axis along its length [8,9]. 

Such a scheme effectively adds a phase acceleration term d to the pendulum equation, 

°° o 

^ = V = 5 -T |a|cos(^ -I- ^), (2.30) 


^The peak-to-peak height of the separatrix is 4y^|a|. 
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Figure 2.7: Tapered PEL phase spaee. An EEL with a linearly tapered undulator has an altered 
pendulum equation and separatrix. If 5/ \a\ = 0, there is no taper and Equation (2.30) reduees 
to Equation (2.16). As 8/\a\ inereases, the separatrix will eontraet, deereasing the area of the 
region in phase spaee with elosed orbits. 


where the phase aeeeleration is given by 



/SK/K h the fraetional rate of ehange of the undulator parameter, and Xs is the dimensionless 
loeation along the undulator where the taper starts. The instantaneous separatrix with a linear 
taper rate has the form [9] 


- Co) + 2 \a\ [sin(C. + <)>)- sin(Co + <)>)], 


where ^0 = 27r — eos ^{—5/\a\) — ^. 

When 5=0 there is no taper, and the separatrix equation reduees to Equation (2.19). As 5 
inereases from zero, the separatrix beeomes distorted as in Eigure 2.7 and the elosed orbit region 
deereases. This reduetion in size of the elosed orbit region ean inerease gain by preventing 
eleetrons from regaining signifieant amounts of energy from the light. When 8/\a\ > 1, there 
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are no closed orbits and electrons cannot be trapped in phase space. Tapering works best when 
the phase acceleration from taper exceeds the maximum untapered deceleration in strong optical 
fields. Estimating the maximum deceleration from the height of the separatrix, we require that 
5 > A^/\a\. When the aforementioned requirements are combined, we see that 

4a/]^ < 5 < |a|. 

Tapering is commonly used in EEL amplifiers, and we will use tapering to optimize an EEE 
design in Chapter 3 [2,9]. 

2.4.4 Short Pulse Effects 

Another practical matter to consider when designing EEEs is the effect of short pulses. Both the 
electron bunch and the light pulses must be synchronized as they travel through the undulator, 
but the electrons are traveling slightly slower than the light. At resonance, the electron bunch 
slips back relative to the light by about one optical wavelength per undulator period. We define 
the slippage distance, NX, as the total distance the electron bunch falls back relative to the 
optical pulse over the length of the undulator. 

If the electron bunch (and the optical pulse) is significantly longer than the slippage distance, 
we do not expect short pulse effects to be significant, because most longitudinal slices of the 
optical pulse will have seen the same charge density over the length of the undulator. On the 
other hand, if the electron bunch length {Op) is on the order of the slippage distance, or 



then we will expect short pulse effects to be significant. 

Eor example, an effect called “optical lethargy” occurs due to uneven amplification of the optical 
pulse over the length of the undulator. Consider the relatively short electron bunch and optical 
pulse in Eigure 2.8. As the electron bunch and optical pulse enter the undulator, the electrons are 
not yet well bunched and are thus not radiating coherently and amplifying the light significantly. 
As the light and electrons travel through the undulator, the electrons begin to bunch, radiate 
coherently, and amplify the light. The strongest amplification of the light occurs near the end 
of the undulator, when the electron bunch has slipped back by . This uneven amplification 
shifts the centroid of the optical pulse, so the light envelope appears to slow down [13,14]. 
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Figure 2.8: Optical lethargy in FEL oscillators. Optical lethargy occurs when the tail of the 
optical pulse is amplified more than the head, causing the centroid of the optieal pulse to travel 
slower than the speed of light. This happens because at the beginning of the undulator (t = 0), 
the eleetrons (shown in red) are not well bunched and are not yet significantly amplifying the 
light (shown in blue). As the electrons start to bunch around T ~ 0.5, they have slipped back 
relative to the optical pulse by half the slippage distance {NX/2) and they amplify the middle 
of the optical pulse moderately. Near the end of the undulator as T 1, the tail of the optieal 
pulse is strongly amplified by the well bunched electrons. After [4]. 


If the electron and optical pulses are exactly synchronized in an FEL oseillator, then over many 
passes the optical lethargy will eause the light pulse to evolve outside of the eleetron pulse. 
Without good overlap with the electrons, the FEL interaction is weak, and there is little to no 
gain. To compensate, we can slightly shorten the optical path length by bringing the resonator 
mirrors eloser together. This technique is referred to as desynchronism and we define the di¬ 
mensionless parameter 

, 2A5 

where A5 is the amount by which the optical cavity is shortened to compensate for optical 
lethargy. Values of A5 are usually quite small—typically less than 5 pm in an infrared FEL 
[14]. The optimum amount of desynchronism varies, but values up to J ~ 0.05 are usually 
viable [9,10]. 

2.4.5 Diffraction 

When we derived the FEL wave equation in Section 2.3.4, we ignored diffraetion in order to 
simplify the derivation. When we include diffraetion, the wave equation gets another term, and 
has the form [9] 

, (2.31) 
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where is the transverse Laplacian and the second term on the right is the same as in Equa¬ 
tion (2.28). Again, this solution assumes that the optical field is slowly varying in amplitude 
and phase along the direction of propagation. 


If we ignore the source term in Equation (2.31), a solution is 


a(v,y, t) = uq 


g-''^/[zo(l+'(r-T„)/zo)] 

1-F/(t-Tvv)/zo 


^0 r-~r^l\zow^(T:)\J<P 
w{x) 


where -\-y^ is the radial distance from the axis squared, w^(t) = 1 -|- (t — XwY Iz^ is the 

square of the optical radius, Xw is the location of the optical waist, zo = ZqIwq is the dimension¬ 
less Rayleigh length,”^ ^ = —im^^[{x — x„)/zQ]+r^{x — Xw)/[z^ + {x — XwY] andao is the field 
amplitude at the optical waist [4]. 


The above solution is the fundamental Gaussian mode; higher-order Gaussian mode solutions 
are also possible [16]. Eow-gain EEE oscillators tend to operate in the fundamental Gaussian 
mode. High gain amplifiers can have significant higher-order mode content. Tilts and shifts of 
the electron beam or cavity mirrors can also lead to optical mode distortion [17]. 


Short Rayleigh Length Oscillators 

In an EEE oscillator, we can control the transverse mode structure using the optical cavity 
mirrors. In particular, we can control the dimensionless Rayleigh length (zq)- When zq is small, 
the beam expands quickly from its waist. This is advantageous for two reasons. Eirst, by 
allowing the optical mode to spread prior to reaching the cavity mirrors, the intensity on the 
mirrors is decreased—an important feature for high power designs [2,18]. Second, since the 
EEE interaction is most pronounced where |a| is large, a single optical wavefront is amplified, 
resulting in good beam quality [2,19]. 

The EEE interaction in the undulator requires good overlap between the electron beam and 
the light. A typical EEE has zq ~ 0.4, which usually allows for the best overlap between the 
light and the electrons over the length of the undulator [19]. The short Rayleigh length (SRE) 
oscillator with zo ~ 0.1-0.2 was proposed in order to reduce the optical intensity on the cavity 
mirrors by allowing the beam cross section to expand more rapidly than longer Rayleigh length 
oscillators. Though the region of overlap between the electron beam and the light is reduced 

^The Rayleigh length Zq is defined as the distance from the mode waist (where the radius is wo by definition) at 
which the cross-sectional area of the beam is doubled. If w(z) is the radius of the mode at z, and Zq is the Rayleigh 

length, then w(z) = w^Jl -^-z} as described in [3,14]. 
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at the ends of the undulator in SRL oseillators, the enhaneed interaetion at the eenter tends to 
dominate [18]. In order to inerease the performanee of SRL oseillators further, the eleetron 
beam ean be foeused in the eenter of the undulator. 
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CHAPTER 3: 

Linear Accelerator Frequency Study 


In this chapter I will describe the study I conducted to determine the effect of accelerator fre¬ 
quency on PEL performance. I begin with a detailed discussion of electron beam generation 
and transport, followed by a motivation for the study. I then describe three models for how 
beam quality is effected by accelerator frequency, the simulations I conducted using these mod¬ 
els, and the results of these simulations. Finally, I consider the implications of this study and 
provide recommendations for future work. 

3.1 Electron Beams 

The goal of electron beam generation for FELs is to create a tightly-packed, relativistic bunch 
of electrons with low energy and angular spreads. The compensation for Coulomb repulsion 
between electrons—and other engineering challenges—makes generation of such quality beams 
more difficult as the bunch charge increases. The quality of an electron beam is determined by 
many factors, and every component in the generation, acceleration, and transport process plays 
a role in shaping the final beam that is delivered to the undulator of an FEE [1]. 

3.1.1 Generation, Acceleration, Transport 

The electron beam begins at a cathode and is rapidly accelerated to relativistic energies by an 
injector into the linac, where it is further accelerated to the required energy of the FEE design. 
Magnetic quadrupoles and dipoles are used throughout the beamline to focus and direct the 
beam. 

Cathode. There are currently two popular cathode technologies used to generate the free elec¬ 
trons that form the electron beams for FEEs: thermionic- and photocathodes. In a thermionic 
cathode, heat provides the energy necessary to eject electrons from the surface of the cath¬ 
ode [2]. The recommended cathode for a USN FEE is a photocathode where incident photons 
from an ultraviolet or visible laser eject electrons by the photoelectric effect [1]. 

Injector. Upon emission from the cathode, the electrons are quickly accelerated to relativistic 
speeds (on the order of 5 MeV) in order to “freeze” their positions relative to one another. This 
freezing effect is a result of relativistic time dilation; in the lab frame, the relativistic electron 
bunch has a slower clock and therefore exhibits slower spreading due to Coulomb repulsion [2]. 
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Accelerator. A linear accelerator is used to accelerate the electron bunches to the ~100MeV 
required for an infrared PEL. In an RF linac, a series of conducting cavities store RF radiation; 
the resulting field gradient accelerates the bunches of electrons as they pass through. 

It is widely agreed that high-power FEFs will require superconducting radio frequency (SRF) 
linacs, as they are more compact, lose less energy in the cavity walls, and can handle higher 
average current. One major disadvantage of SRF linacs is that they require a cryoplant and 
must operate at temperatures as low as 2 K [2]. 

Transport, focusing. Dipole and quadrupole magnets are used to direct and focus the elec¬ 
tron beam throughout the beamlines of FEFs. For example, prior to entering the undulator 
of a short Rayleigh length (SRF) oscillator, the electron beam is focused to the center of the 
undulator—using quadrupole magnets—where the optical mode is also focused, enhancing the 
FEE interaction there [18-20]. In an FEE amplifier, the electron beam is typically “matched” 
so that its average radius remains constant as it passes through the undulator; this is necessary 
due to the longer undulators required by amplifiers. Quadrupole focusing will be discussed in 
detail in Chapter 4. 

3.1.2 Characteristics 

From the perspective of FEE performance, an electron beam has many important characteristics 
that can be quantified. 

Energy. The electron energy is typically specified either directly in electron volts or by the rel¬ 
ativistic Forentz factor 7 , as discussed in Chapter 2. An FEE operating at infrared wavelengths 
will require an electron beam of ~100MeV with a corresponding Forentz factor of 7 ~ 200. 

Frequency. In RF linacs, electron beams are not continuous streams of electrons—they are 
pulsed so that the electrons are delivered in bunches. The pulse repetition frequency (prf) must 
either be equal to or a submultiple of the accelerator frequency in the cavities in order for the 
bunches to experience the appropriate RF phase for acceleration [21]. 

Bunch charge. The average current in an electron beam is determined by the amount of charge 
in each bunch (pulse) and the prf. For the same average current, an electron beam with a lower 
prf will have more charge per bunch, and vice versa. High-power FEFs require bunch charges 
on the order of a nanocoulomb. 
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Pulse length (duration). Each bunch of electrons is spread over some distance along the direc¬ 
tion of travel; this spread is referred to as the pulse length or duration. Together with the bunch 
charge, the pulse duration (and shape) determines the peak current of the electron beam. The 
pulse duration for a high-power EEL is usually on the order of a picosecond. 

Radius. The radius of the electron beam is measured perpendicular to the direction of travel, 
and is usually about a tenth of a millimeter. Smaller radii increase the Coulomb repulsion 
forces between the electrons but can also strengthen the EEE interaction, resulting in increased 
performance. 

Emittance. Beam quality is quantified by emittance. Eongitudinal emittance is the product of 
the bunch length and energy spread and is typically specified in units of keV-ps. Transverse 
emittance is the product of the beam radius and the angular spread in the electron velocities; it 
is typically given in units of mm-mrad [8,21]. 

3.2 Motivation 

In the preferred electron beam source for naval application—the SRE linac—electrons are ar¬ 
ranged in bunches with lengths on the order of picoseconds which are typically separated from 
each other by a few nanoseconds [1]. As discussed above, the accelerator frequency is the 
frequency at which the electric fields in the superconducting cavities vary, and this frequency 
determines the maximum repetition rate at which the electron bunches can be delivered. 

There are engineering trade-offs associated with any choice of accelerator frequency. At present, 
high-power EEE designs utilize frequencies from a few hundred megahertz to well over one gi¬ 
gahertz [5]. Given an average current requirement, higher frequencies have lower bunch charge 
are require smaller accelerating cavities. Eower frequencies result in more charge per bunch 
and require less refrigeration; in fact, lower frequency designs might be operated at 4 K (higher 
frequency SRE linacs require refrigeration to 2 K) [22]. 

We will not consider the benefits of a higher or lower accelerator frequency from an engineering 
standpoint in great detail—that is beyond the scope of this thesis. It is reasonable to postulate, 
however, that designing an EEE system with an SRE linac that operates at 4 K instead of 2 K 
(twice the absolute temperature) could reduce the overall size, cost, complexity, and power 
requirements of the EEE system [2]. As such, it would be prudent to consider the effects of 
accelerator frequency on EEE performance to reveal any trade-offs and inform future design 
decisions. 
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3.3 Emittance Scaling Models 


To perform a meaningful comparison between the different accelerator frequencies that we wish 
to investigate, we must hold as many variables constant as possible. Most importantly, we will 
need to keep the peak and average currents fixed, as current plays a critical role in determining 
PEL performance. As we change the frequency, this will require us to adjust the bunch charge 
accordingly. Doing so will have some effect on the beam emittance, so we will need some sort 
of model for how the emittance will scale with accelerator frequency. 

Scaling laws for emittances are difficult to derive analytically (or determine numerically). With 
the assistance of Prof. T. Smith at Stanford University, we have developed three simple models 
of how the electron beam’s emittance might scale as the frequency of the accelerator is changed, 
given a fixed average and peak current. 

In Figure 3.1 we outline the rationalization behind these models. Given an electron bunch of 
charge Q from an operating FEL/linac pair with bunch length L, radius R, accelerator frequency 
/, transverse beam emittance e± = and longitudinal beam emittance £/ = (Figure 3.1a), 
we note that if the accelerator frequency is halved, we must double the charge per bunch in order 
to keep the average current constant. We can lengthen the pulse and keep the charge density 
the same (which satisfies the constraint that we keep the peak current fixed) by placing two 
bunches, each with charge 2, end to end. This will keep the average current fixed with a new 
accelerator frequency //2. The transverse beam emittance should be unaffected by this change, 
so we must determine how the longitudinal emittance will be affected. 

In the worst case, we might expect the longitudinal emittance to double if the bunch length dou¬ 
bles and the average energy spread remains the same, as shown in Figure 3.1b (this corresponds 
to our model where £/ 1//). However, it is reasonable to assume that the electrons in the 

middle of the combined bunch 2Q will be distributed in energy space as shown in Figure 3.1c, 
resulting in a longitudinal emittance £/ where £/ < £/ < 2£;. 

The rationalization in Figure 3.1 motivated us to run our simulations using three models for the 
scaling of longitudinal emittance with accelerator frequency. In all of the models we assume 
that the transverse emittance is unaffected by changing the accelerator frequency. In the first 
model—the “worst case” model—we assume that the longitudinal emittance scales inversely 
with the frequency. This model corresponds to Figure 3.1b. In the second model—the “com¬ 
promise” model—we assume that the longitudinal emittance scales inversely with the square 
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Figure 3.1: Rationalization of the longitudinal emittance scaling models. The longitudinal 
emittance of a beam is related to the energy spread of its constituent electrons and will vary 
with the accelerator frequency. Here we outline the rationale behind our longitudinal emittance 
scaling models by showing how the energy spread of an electron bunch might change when the 
accelerator frequency is halved. After [22]. 


root of the frequency. This model corresponds to the case shown in Figure 3.1c. We refer to 
this model as the “compromise” model because in the third model—the “best case” model—we 
assume that the longitudinal emittance of the beam is independent of the accelerator frequency. 

The worst case and best case models provide a lower and upper bound estimate on the per¬ 
formance effects we should expect from changing the accelerator frequency. The compromise 
model is what we consider to be the most realistic case, and is based on a physical model 
developed by Prof. T. Smith [22]. 


3.4 Amplifier Simulations 

A previous study conducted by Prof. J. Blau, Prof. W. B. Colson, Prof. K. Cohn, Prof. T. Smith, 
and Matt Stanton investigated the effect of varying the accelerator frequency on FEL amplifiers. 
In that study, amplifier simulations were run at three beam energies using the compromise model 
we described above. I extended that study to include the best case and worst case models. 
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£/ (keV-ps) 

/ (MHz) 

1 // 

1 / ^/f constant 

300 

320 

226 

160 

400 

240 

196 

160 

500 

192 

175 

160 

600 

160 

160 

160 

700 

137 

148 

160 

800 

120 

139 

160 

900 

107 

131 

160 


Table 3.1: The longitudinal emittances used in the oscillator and amplifier simulations. A ref¬ 
erence beam quality of 160 keV-ps was chosen at 600 MHz, and the emittances at the remaining 
frequencies were calculated using each of the three models described in Section 3.3. The results 
of the simulations are summarized in Figure 3.3 and Figure 3.5. 


Simulations were run for frequencies from 300 MHz to 900 MHz with an electron beam energy 
of 100 MeV, average beam current of 80 mA, and a peak current of 1 kA. A reference electron 
beam with a longitudinal emittance of 160keV-ps was chosen at 600 MHz, the middle of the 
RF range considered in the study. The emittances at lower and higher frequencies were then 
computed using each of the models described in Section 3.3. The longitudinal emittances that 
were used in the simulations are listed in Table 3.1. The undulator taper rate, taper start location, 
and initial electron phase velocity were all optimized for each accelerator frequency. 


3.4.1 4D Amplifier Simulation Programs 

Prof. J. Blau and Prof. W. B. Colson have developed a suite of FEL simulation software at NPS. 
One of the programs, f el4dl (and its optimization variant, f el4dlv) provides 4D simulation of 
an FEL amplifier; it simulates the optical wavefronts and electrons in a fully three-dimensional 
{x, y, z) space as they evolve in time. The f el4dl and f el4dlv programs were used to run the 
FEL amplifier simulations for this study. 

Both f el4dl and f el4dlv track the evolution of the optical wavefronts using the electromag¬ 
netic wave equation (derived from Maxwell’s equations), and they follow the evolution of the 
electrons using the Lorentz force equation as described in Chapter 2. 

The programs use parallel computing to reduce run time, executing on ~100 CPU cores at a 
time. In f el4dl and fel4dlv, each core handles the simulation of a longitudinal “slice” of 
the optical pulse. Each slice is stored on a transverse grid of around 200 x 200 points and is 
populated with approximately 30,000 sample electrons. 
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The transverse grid and longitudinal sliees allow the programs to model sophistieated effects 
such as transverse optical modes and slippage. For example, as each slice of the optical pulse 
travels down the undulator, the slower-moving electrons are successively handed off to the 
following optical slice. By the end of the undulator, the electron bunch will have slipped back 
by one slippage length (NX) relative to the optical pulse. 

The simulations require a plain text input file with the desired parameters of the FEL design be¬ 
ing tested. Upon completion, the simulation is summarized in a graphical format. An example 
of the program output is shown in Figure 3.2. The output includes a list of dimensionless phys¬ 
ical parameters and multiple plots which characterize the performance of the FEL. In addition 
to the graphical output that is shown in Eigure 3.2, text output is provided, including: numerical 
parameters used in the simulation, final gain, extraction, energy spread, and other diagnostic 
values. 

The f el4dlv program is an optimization variant of f el4dl which uses a nearly identical in¬ 
put file, with one exception: the f el4dlv input file includes a parameter that indicates which 
variables should be optimized. This allows both the taper rate and taper start location in the 
undulator to be optimized with very little user effort. Other variables that can be optimized 
include the initial electron phase velocity and the spread in phase velocities due to the electron 
energy spread. 

The top section of the output shown in Eigure 3.2 includes the dimensionless physical parame¬ 
ters used to run the simulation. In the plots, a color bar ranging from dark blue to cyan is used 
to indicate, depending on the plot, amplitude or value: dark blue is used for an amplitude or 
value of zero and cyan is used for a maximum amplitude or value. The top section of the output 
includes four column plots; each column has a horizontal axis corresponding to one of x, y, z, 
or V and shows the evolution of the electron beam and optical mode for T = 0 —1. The top 
of each column contains a plot at T = 1, while the bottom has a plot at t = 0. The first two 
columns show cross section slices in x and y of the optical amplitude (as a color plot) and the 
electron beam envelope (in red). The third column shows the optical and electron pulse shapes 
in the longitudinal (z) direction, and the final column shows a power spectrum (in v). 

The bottom row of graphs in the output begins on the left with an electron phase velocity 
distribution plot for T = 0 —)■ 1. In the next plot, the final electron phase space (at T = 1) is 

O 

plotted in the dimensionless ^ and V = C coordinates. Often it is easy to observe electron 
bunching and energy loss (or the lack thereof) here. The last two graphs show a surface plot 
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-Jt/2 C 371/2 -3x3 -12 X 12 


Figure 3.2: Sample 4D simulation program output (f el4dl). This output was generated for the 
600 MHz oseillator with an assumed £; 1 /^/J dependeney. A list of dimensionless physieal 
parameters is provided at the top. In the middle, the electron bunch and optical pulse structure 
are shown in x, y, and z along with a power spectrum (in v). At the bottom, an electron phase 
velocity distribution, phase space plot, and optical mode surface plots are shown. 


of the transverse optical amplitude at the end of the undulator and at the first optical element 
beyond the undulator. 

The f el4dl output provides all the information required to evaluate the performance of a given 
PEL amplifier configuration. In particular, we are concerned with extraction of light and final 
energy spread in the electron beam. Higher extraction allows for a lower average current in the 
electron beam for a desired output power and a lower energy spread allows for more efficient 
energy recovery from the electron beam [2]. 
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3.4.2 Computing 

The amplifier portion of this study required approximately 1800 eomplete simulations. At eaeh 
aeeelerator frequeney, many simulations were run to optimize the initial eleetron phase veloeity, 
undulator taper rate and taper start loeation. Eaeh simulation required ~250 minutes of CPU 
time on a eluster eomputer. 

We ran the amplifier simulations on the EEL group eluster, phf el. ern. nps . edu (in the original 
study, the simulations were run on the NPS eluster, hamming.uc .nps . edu). The simulations 
were run over a period of several weeks, and required over 7500 CPU hours to eomplete. 

3.4.3 Amplifier Results 

The results of the amplifier simulations are plotted in Eigure 3.3. The red line (eireles at eaeh 
data point) represents the extraetion versus aeeelerator frequeney when the longitudinal emit- 
tanee seales inversely with the aeeelerator frequeney; this was the worst ease model. Even in 
the worst ease, we noted that the amplifier performanee inereased slightly as we lowered the 
aeeelerator frequeney. The performanee inereases are attributable to redueed energy spread and 
pulse slippage effeets (a benefit of the longer pulse) in all three eases. 

The green line (triangles) represents the amplifier extraetion versus aeeelerator frequeney when 
the longitudinal emittanee seales inversely with the square root of the aeeelerator frequeney. 
This model is the most realistie and suggests that signifieant performanee gains eould be real¬ 
ized by redueing the aeeelerator frequeney of a given system. Over the range of frequeneies 
ineluded in this study, the extraetion appears to inerease linearly as the aeeelerator frequeney 
is lowered. The blue line (squares) represents the extraetion versus aeeelerator frequeney when 
the longitudinal emittanee is held eonstant and represents an upper bound on the performanee 
inerease that eould be realized by redueing the aeeelerator frequeney of a given EEE. 

3.5 Oscillator Simulations 

Simulations were also run for an EEE oseillator to determine the effeet of ehanging the aeeeler¬ 
ator frequeney on performanee. As in the amplifier study, simulations were run for frequeneies 
from 300 MHz to 900 MHz with a beam energy of lOOMeV, average beam eurrent of 80 mA, 
and a peak eurrent of 1 kA. The same referenee eleetron beam that was used for the amplifier 
study was used for the oseillator study—a beam with a longitudinal emittanee of 160keV-ps 
at 600 MHz. The emittanees at lower and higher frequeneies were then eomputed using eaeh 
of the models deseribed in Seetion 3.3 and are listed in Table 3.1. The beam radius, Rayleigh 
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Figure 3.3: Predicted PEL amplifier extraction versus accelerator frequency for the three longi¬ 
tudinal emittance models. The extraction increases as the accelerator frequency decreases in all 
three models. The performance gains are largely due to reduced energy spread and slippage ef¬ 
fects as the electron bunch length increases. All three models have identical performance at the 
600 MHz reference frequency since the electron beam emittance was set to 160keV-ps there. 

length, optical cavity quality factor, number of undulator periods, and optical intensity on the 
cavity mirrors^ were all optimized for each accelerator frequency. 

3.5.1 3D Oscillator Simulation Program 

The wavevnm program provides 3D simulation of an EEL oscillator; it simulates the optical 
wavefronts and the electron bunch in the transverse (x,y) dimensions and time. The wavevnm 
program was used to run the EEE oscillator simulations for this study. Eike fel4dl and 
fel4dlv, wavevnm follows the evolution of the optical wavefronts and the electrons using 
Maxwell’s wave equation and the Eorentz force equation. 

Unlike the 4D programs, wavevnm follows only a single longitudinal slice of the optical pulse 

*An engineering limitation requires that the optical intensity on the mirrors be less than ~100kWcm^^ to 
prevent damage or distortion to the mirrors [1,2]. 
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and electron bunch. This amounts to an assumption that most slices of the optical pulse are iden¬ 
tical and see roughly the same current density. This assumption is valid for cases where the pulse 
length is long compared to the slippage distance, as it is for the PEL oscillators studied here. 

The simulations use parallel computing to reduce run time, executing on many CPU cores at a 
time. In wavevnm, each core runs the simulation with a different average initial electron phase 
velocity, Vq. The electron phase velocity corresponds to an optical wavelength, so running the 
simulation for multiple values of Vq is equivalent to running the simulation with multiple optical 
wavelengths. In an actual EEL, mode competition determines which value of V (which optical 
wavelength) dominates. The simulation requires that the user examine the output to determine 
which value of Vq resulted in the greatest extraction; this would be the dominant wavelength in 
an actual EEE.^ 

The simulation requires a plain text input file with the desired parameters of the EEE design 
being tested. Upon completion, the simulation is summarized in a graphical format. An example 
output is shown in Eigure 3.4. The output includes a list of dimensionless physical parameters 
and multiple plots which characterize the EEE performance. In addition to the graphical output 
that is shown in Eigure 3.4, text output is provided, including: numerical parameters used in the 
simulation, final gain, extraction, energy spread, and other diagnostic values. 

The upper right section of the output shown in Eigure 3.4 includes the dimensionless physical 
parameters used to run the simulation. In the plots, a color bar ranging from dark blue to cyan is 
used to indicate, depending on the plot, amplitude or value: dark blue represents an amplitude or 
value of zero and cyan represents a maximum amplitude or value. In the upper left, a cross sec¬ 
tion of the optical wavefront at the output mirror is shown versus the pass number. To the right 
of this, the optical wavefront is shown as a surface plot at the output mirror on the final pass. 

In the second row of plots, a cross section of electron positions, optical amplitude (in the un- 
dulator) and optical mode shapes (at the mirrors) are all plotted on the final pass. A random 
sampling of electrons is shown as red dots and mode shapes are plotted on both ends of the 
undulator (a fundamental Gaussian is shown in magenta and the optical mode shape for the 
simulation is shown in yellow). 


4D oscillator simulation is able to run just one simulation which self-consistently evolves to the dominant 
wavelength. It is worth noting that in an amplifier, it is not necessary to consider very many initial values of Vq, as 
the seed laser is the dominant factor in determining the output wavelength. 
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FEL 3D SIMULATION 
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Figure 3.4: Sample 3D simulation output (600 MHz oscillator, £/ oc A list of dimen¬ 

sionless physical parameters in provided in the upper right. In the upper left, a plot of the 
optical mode amplitude at the output mirror is shown. In the middle, the electron beam and 
optical mode are shown in the undulator. At the bottom, an electron phase velocity distribution, 
phase space plot, power in the optical cavity, and gain plots are shown. 


The last row of plots begins with an electron phase velocity distribution versus the pass number. 
The next plot is a phase space plot of a random sampling of the electrons in the dimensionless 

O 

^ and V = ^ coordinates. In the lower right, the optical power in the cavity and optical gain are 
both plotted versus pass number. 

The wavevnm output provides all the information required to evaluate the performance of a given 
FEL oscillator configuration. In particular, we are concerned with extraction and mode quality. 
Higher extraction allows for a lower current in the electron beam for a desired output power and 
mode quality affects the propagation of the laser beam through the atmosphere to the target. 
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3.5.2 Computing 

The oscillator portion of this study required approximately 240 complete simulations. At each 
accelerator frequency, many simulations were run to manually^ optimize the resonator cavity 
output coupling factor, number of undulator periods, optical Rayleigh length, initial electron 
phase velocity, and electron beam waist radius. Each simulation required ~240 minutes of 
CPU time on a cluster computer. The oscillator simulations were run on the EEL group cluster, 
phf el. ern. nps . edu. The simulations were run over a period of several weeks and required 
around 960 CPU hours to complete.^ 

3.5.3 Oscillator Results 

The results of the oscillator simulations are plotted in Eigure 3.5. As with the amplifier re¬ 
sults, the red line (circles at each data point) represents the extraction versus accelerator fre¬ 
quency when the longitudinal emittance scales inversely with the accelerator frequency, the 
green line (triangles) represents the extraction versus accelerator frequency when the longitudi¬ 
nal emittance scales inversely with the square root of the accelerator frequency, and the blue line 
(squares) represents the extraction versus accelerator frequency when the longitudinal emittance 
is held constant. 

It is immediately apparent from Eigure 3.5 that EEE oscillators are not as sensitive to changes 
in the accelerator frequency as amplifiers. Even in the best case model (in blue) the extraction 
increases only slightly as we reduce the accelerator frequency. EEE oscillators are not as sensi¬ 
tive to beam quality as amplifiers due to their shorter undulators. Though the increase in pulse 
length as the accelerator frequency is lowered does reduce the negative effects of energy spread, 
the effect on the extraction is not as marked as it is for amplifiers. 

3.6 Conclusion 

Motivated by the potential engineering advantages of lowering the frequency of an EEE linac, I 
conducted a study to determine the effect this would have on EEE performance. I used several 
models for how longitudinal emittance scales with accelerator frequency and applied each of 

^At the time of this writing, there was no automatic optimization variant of wavevnm available. 

‘^The astute reader might notice a seeming discrepancy in the CPU time and number of simulations required 
to complete the oscillator simulations versus the amplifier simulations. In the oscillator study, most of the opti¬ 
mization work was done “by hand,” reducing the number of simulations required, whereas the amplifier program 
f el4dlv optimizes by brute force, running many more simulations than is strictly necessary (but requiring less 
offline calculation). 
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Figure 3.5: Predicted PEL oscillator extraction versus accelerator frequency for the three lon¬ 
gitudinal emittance models. In the worst case (red line), the extraction is constant; in the other 
cases, the extraction increases slightly as the accelerator frequency decreases. All three mod¬ 
els have identical performance at the 600 MHz reference frequency since the electron beam 
emittance was set to 160keV-ps there. 

those models to simulate high-current EEL designs in both oscillator and amplifier configura¬ 
tions using the NPS FEE simulation suite. 

The amplifier simulations suggest that significant performance gains might be realized by using 
a lower accelerator frequency. In what is the most realistic emittance model, EEE performance 
increased linearly as the accelerator frequency was lowered from 900 MHz to 300 MHz. The 
performance gains were driven largely by a lengthening electron pulse and the commensurate 
reduction in energy spread and slippage effects. 

Oscillators are not as sensitive to electron beam quality as amplifiers, so it was not surprising 
that the performance effects of changing the accelerator frequency were not as dramatic in the 
oscillator simulations. The simulations were 3D, so they did not account for short pulse effects, 
but short pulse effects are not as important in oscillators as they are in amplifiers. The oscillator 
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designs performed well at all of the frequencies in this study, so the engineering advantages 
of a lower accelerator frequency alone should provide motivation to consider building PEL 
oscillators that use lower accelerator frequencies. 

In summary, I believe that lowering the accelerator frequency of both EEL oscillators and am¬ 
plifiers (particularly those utilizing SRE linacs) is worthy of consideration due to the potential 
for improved performance and reduction in size, complexity, cost, and power requirements of 
the linac system. Both the oscillator and amplifier designs that I simulated had more than suffi¬ 
cient performance at lower frequencies to justify the design of new EEEs that utilize SRE linacs 
with lower frequencies. 

3.7 Future work 

There are opportunities for additional work to be done to extend this study, such as running 
the simulations for amplifiers and oscillators at other wavelengths (such as optical, terahertz, 
and x-ray) to verify the conclusions of this study at those wavelengths. Running the oscillator 
simulations on the 4D oscillator code that was recently developed by Prof. J. Blau at NPS could 
provide validation to both the 4D oscillator code and these 3D results. It could also be beneficial 
to run start-to-end simulations using energy recovery linacs and optimizing the electron beam 
transport using codes like PARMEEA, EEEGANT, GPT, DIMAD, and EEESIM. With these 
start to end simulations, we could more accurately determine how the accelerator frequency 
affects the longitudinal and transverse emittance. 
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CHAPTER 4: 

Electron Trajectories in Quadrupoles and an Undulator 


Every electron beam will have a spread in electron positions and velocities transverse to the 
axis of the beamline. Without external focusing, the electron beam will inevitably expand ra¬ 
dially due to the spread in transverse velocities, decreasing the electron density and potentially 
resulting in electrons striking beamline components. 

This is undesirable for many reasons, including reduced FEE gain, potential damage to beam¬ 
line components, and the production of unwanted x-ray and background radiation [2]. Magnetic 
quadrupoles, like those shown in Figure 4.1, are used to focus the electron beam and prevent 
spreading as the bunches make their way from the cathode through the various beamline com¬ 
ponents and to the FEE undulator [21]. 



Figure 4.1: A quadrupole triplet in the beamline of the x-ray FEE at SACEA in Japan. 
Quadrupole focusing is essential to keep the highly energetic electrons confined to the beam 
pipe on their long journey to the undulator section. The yellow structure shown here supports 
each quadrupole, and the brown coils are the electromagnets that produce the quadrupole field. 
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In the undulator of an FEL, the eleetrons undergo the betatron oseillations introdueed in See- 
tion 2.4.2. These oseillations eause a different type of foeusing, and I will deseribe them in 
detail in Section 4.2. 

Section 4.4 describes the simulation that I developed to track the trajectories of electrons as they 
are focused by magnetic quadrupoles and the FEL undulator. I will conclude this chapter by 
discussing potential applications of the simulation and opportunities for future work. 


4.1 Quadrupoles 

The fields of a magnetic quadrupole are created with four hyperbolic pole faces arranged as 
shown in Figure 4.2. There are two planes of symmetry created by this arrangement, and the 
approximate fields near the quadrupole axis are given by 

B, = Bq^- (4.1a) 

a 

B,=Bq- (4.1b) 

a 

where a is the radial distance from the axis to the pole face as shown in Figure 4.2 and Bq is the 
component of the field strength measured when x = a (or y = a) [2\]. 



Figure 4.2: Fields of a magnetic quadrupole. The fields have two planes of symmetry: the 
y = X plane and the y — —x plane. An electron traveling in the z direction (into the page) will 
be focused in the y direction and defocused in the x direction. After [21]. 
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It is also possible to use axially symmetric electrostatic or magnetic focusing fields, such as 
those produced by an electrostatic circular aperture or a magnetic solenoid, but the focusing 
strength of their fields is usually an order of magnitude weaker for relativistic beams than the 
focusing strength produced with magnetic quadrupoles [ 21 ]. 


4.1.1 Equations of Motion 

To derive the equations of motion for electrons in a quadrupole with the fields of Equation (4.1), 
we begin with the relativistic Lorentz force equation: 


d(rfi) 

dt 


^(E + J3 xB 

me 


In a magnetic quadrupole the electric field is zero, and magnetic fields do no work, so 7 = 0 and 
we have (with 5^ = 0 in an ideal magnetic quadrupole) 






ymc 

d 

ymc 


/5xB, 

+ ftBxy + (ftBy - ftBj) 2]. 


In an EEL, q = —e and the electrons are highly relativistic, so we use the approximations ~ 1 
and Py ^ /3z and we get 


Px ~ 

- By 

ymc 

Py- 

1 

ymc 

V'x ~ 

— By 

ym 


——B 

ym 

X ~ 


y ~ 

cBq^ 


aym 

aym" 

X ^ 

(OqX 

y ~ 

-(oly 


(4.2) 


where COq = eBQ/aym and we substituted the quadrupole fields from Equations (4.1a) and 
(4.1b). The possible solutions arex = Acosh((Ogt) + 5 sinh(( 0 gt) andy = Ceos(tOgt) +Dsin((Ogt) 
where A, 5, C, and D are constants that will depend on initial conditions. Thus, this quadrupole 
is focusing in the y direction and defocusing in the x direction [ 21 ]. 


4.1.2 Quadrupoles as Lenses 

It is often possible to treat the quadrupole as a thin lens that acts more or less instantaneously. 
Before showing this, we convert Equations (4.2) to dimensionless time; define l/dt^ = c^jL^dt^ 
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and write 


Sdq 

X 

ajm 

oo 

X ^ 

aymc^ 

OO 

X KX 


y 


A? 

y rs-/ 


eBQ 

— 

aym 

aymc^ 


OO 

y ~ -Ky, 


where Lq is the width of the quadrupole field and k = bLqBq! aymc^ . 


(4.3) 


To show how the quadrupole can act like a lens, we consider an electron with a transverse 
velocity (which we have assumed to be very small compared to its longitudinal velocity). As 
the electron enters the quadrupole, it will obey Equations (4.3), so we can write 


/ ° 7 

Vxf ~ Vxi + y VxdT 

oo 

^Vxi + / xdx 

.h, 

C'^e 

^Vxi + / KxdX, 
Jts 


where Vxj and Vxi are the velocity before and after entering the quadrupole, respectively, and 
Xs and Xe are the dimensionless times corresponding to the start and end of the quadrupole 
fields.^ We have assumed that Vx is small and the width of the quadrupole field is small, so x 
is more or less constant for the duration of the transit through the quadrupole fields. With this 
approximation we get 


Vxf ~ 


Vxi + / Kxdx 

./Tj 

Vxi T K(^Xe Xs)x 
Vx/ + (kAx)x. 


where Ax = Xg — Xs- We define II f = kAx and write 


Vxf ~ Vxi + y 


Vyf 


Ai 


/■ 


(4.4) 


'clearly the fields do not start and end abruptly as we imply here. However, for a first order approximation, 
this calculation will suffice. For any given field, a hard-edge approximation can be found [21]. 
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Equations (4.4) describe the action of a lens with focal length / that is focusing in y and defo- 
cusing in x (where / > 0 in this case since kAt is always positive). 


4.1.3 Quadrupole Doublets and Triplets 

Two quadrupoles with focal lengths /i and /2 separated by short drift space of length s form a 
quadrupole doublet with an equivalent focal length given by 


1 _ 1 1 5 


(4.5) 


If the quadrupoles are of equal strength and oriented so /i = / and /2 = —/, then Equation (4.5) 
becomes 


1 _ 5 

F~p' 


(4.6) 


A quadrupole triplet is made up of a quadrupole with a focal length / surrounded by two 
quadrupoles, each with a focal length of —2/ [21]. 


4.2 Betatron Motion 

As discussed briefly in Chapter 2, in addition to the fast “wiggling” motion caused by the 
transverse periodic magnetic field in the undulator of an EEL, the electrons also undergo a 
slower betatron motion which is caused by the transverse variation in the magnetic field strength 
of the undulator. In a linear undulator, for example, the field strength will increase as the 
electron gets closer to the magnetic pole faces. Consider the field of such a linear undulator: 

B = 5osin(fcoz)cosh(fcoy)y + 5ocos(fcoz)sinh(fcoy)z, (4.7) 

where the undulator is oriented so the pole faces are in the y direction. If we expand Equa¬ 
tion (2.1) assuming that there is no light in the undulator (so y = 0), we get 

= —^ [{^yBz - My) X - Mzf + My^] , (4.8) 
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where By and are the components of Equation (4.7). We get /3x by integrating the x- 
component with respect to time: 


/3x = 


-sBq 

ymc^ 


ycos{kQz) sinh(fcoy) -zsin(fcoz) cosh(fcoj)], 


—cBq 1 d 
ymc^ ko dt 


cos(fcoz) cosh(fcoy)], 


Px = T~^ cos(koz) cosh(koy) + C, 
koymc^ 


where C is a constant of integration. For “perfect” injection into a sinusoidal path, C = 0, and 
we will assume that this is the case for simplicity. We will also assume that the excursions made 
by the electrons from the axis of the undulator are small so kgy 1. The undulator parameter 
for a linear undulator is 

^ (^BfffigXQ £Bq 

iTtmc^ y/2komc^ 

where B^^s = Bo /for a linear undulator, Bq is the peak magnetic field strength, and the 
undulator period is Aq = 2n/ko. Using K, we simplify the expression for Px to 


I5x ^ — ^^ cos(koz) ~ — ^^ cos(koct), 
where we made the approximation that z ct. We integrate again to get .r: 


(4.9) 


Px^ 



V2K 

r 

VlK 

r 

V2K 

hr 


cos{koct) 

r' 

-cos{koct) 
sin(koct) +.ro- 


Where xq is the x displacement at t = 0. We assumed perfect injection into a sinusoidal path, so 
xq = 0. Recall that ko = 2 k/ h, so we can write 


X ^ 


KXq 

\/2Ky 


sin(A:oct). 


(4.10) 


This is not the slow betatron motion (the argument of the sine function is koct )—this is the fast 
“wiggling” motion that produces light. We will need this expression to solve for the betatron 
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motion, however. The y-component of Equation (4.8) is 


= -fco4cos(fcoz) sinh(fcoy) 

7 

—2cK^ 

— -2—A:oCOS^(A:oz)sinh(A:oy)cosh(A:oy) 

cos^(A:oz) sinh( 2 A:oy). 


r 

—cK^ko 


T 


The average value of cos^{koz) over a period is 1/2, so the average value of the previous ex¬ 
pression is 

Again, we assume that the excursions made by the electrons from the undulator axis are small 
(koy 1 ) so sinh( 2 A:oy) ~ 2 A:oy and we can write 



c^K^kl 


which is a differential equation describing simple harmonic motion. We convert to dimension¬ 
less time (t = ct/L) and get 


OO 

7 = - 



(4.11) 


where we have defined the dimensionless betatron frequency as 


®/3 


Kk^L _ IkNK 

y ~ y 


(4.12) 


In a typical EEL, N ~ 100, ^ ^ 1, and y ~ 100, so betatron frequency values are on the order 
of 2k. This corresponds to one oscillation over the length of the undulator (t = 0 —)■ 1). Solve 
for y(T) by integrating twice with respect to dimensionless time, and we get 


7('r) 


yocos{cOpT) -h 


L0y . 

—^ sm 

(Op 


{(Opx), 


(4.13) 


where yo is the displacement at T = 0, By is the angular displacement of the electron velocity 
from the undulator axis at T = 0, and LOy is the rate of change of y with respect to dimensionless 
time at T = 0 [10,15]. We can get a feeling for the magnitude of these betatron oscillations 
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compared to the fast wiggling motion by inserting some typical values into Equation (4.10) and 
Equation (4.13). If we take ^ = 1, = 100, Aq = 2cm, y = 100, yo = 1 and estimate the 

angular spread as Oy ^ 'Ijq/L, we get the maximum amplitude of the fast wiggling motion to be 


-^max 


(l)(2cm) 

(V2;r)(100) 


50 pm. 


To estimate jmax, the maximum amplitude of the slow betatron motion, we note that the factor 
multiplying the sine in Equation (4.13), 

dy 2yo yoT Jo 

is about three times smaller than yo (which multiplies the cosine in Equation (4.13)), so we 
approximate ymax ~ 1 mm. It appears that the betatron oscillations are roughly an order of mag¬ 
nitude greater in amplitude than the fast wiggling motion. 


4.3 Previous Work 

The simulation I developed is based on two existing programs: quad and betaxy. They were 
written by Prof. W. B. Colson and Prof. J. Blau using the C programming language, use plain 
text input files or interactive keyboard input, and generate an output graphic that summarizes 
the results of the simulation. 

4.3.1 Quadrupole Simulation 

The quad program simulates electron trajectories through a single quadrupole that is modeled 
as a thin lens. The program requires two sets of input parameters: one to describe the electron 
distribution and another to describe the quadrupole. The initial electron distribution is assumed 
to be Gaussian with zero mean in both position and angular spread in the x and y directions, and 
has no longitudinal structure (all electrons are assumed to lie in the same x-y plane). There are 
thus four input parameters for the electron beam: 

• Ox, the standard deviation of the electron x-coordinates; 

• Oy, the standard deviation of the electron y-coordinates; 

• Oq^, the standard deviation in the angular spread in the x direction; and 

• O 0 ^, the standard deviation in the angular spread in the y direction. 
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These parameters are speeified at T = 0 (the starting point of the beam in the simulation). All 
input parameters are normalized to be dimensionless. 


1 


x(x) 


-1 

1 


y (X) 

-1 

Figure 4.3: Sample output of the quadrupole simulation program (quad). A single quadrupole 
with a foeal length / = 0.5 is plaeed at Tf = 0.7. The eleetron distribution has an angular 
spread of — Oq^, = 0.05 and a transverse spread of (7x = CJ); = 0. 1 . The eleetron trajeetories 
are plotted in both jc and j as T = 0 —)■ 1. 

The quadrupole focusing strength is specified by entering the focal length, /. A positive value 
for / produces focusing in x and defocusing in y, and vice versa when / is negative. The position 
is indicated with Xf, the dimensionless time at which the electrons encounter the quadrupole. 

The program produces a graphical output containing the results of the simulation, an example 
of which is shown in Figure 4.3. The output includes a list of the physical input parameters 
as well as plots of the electron trajectories x{x) and y(T). Note that the quadrupole placed at 
Xf = 0.7 with a focal length of / = 0.5 is focusing in x and defocusing in y. 

4.3.2 Betatron Motion Simulation 

The betaxy program simulates the betatron oscillations of electrons as they pass through an 
FEL undulator. Like quad, the program requires two sets of input parameters: in betaxy, one 
describes the electron distribution and another describes the undulator. The electron distribution 
is described using the same input parameters as quad (along with the same assumptions regard¬ 
ing its electron distribution). An additional parameter, xp, indicates the position in the undulator 
at which the electron beam distribution parameters are given. The undulator is described by its 


*** Quadrupole Focusing *** 

f=0.5 a^=0.1 ®ex"°• 

x_=0.7 a =0.1 Gr, =0.05 

f y by 



0x1 
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betatron frequency, which was defined in Equation (4.12). Both and can be specified, 
so there can be unequal betatron motion in x and y. 


y 


Figure 4.4: Sample output of the betatron motion simulation program (betaxy). The undulator 
in this simulation has a betatron frequency of = n. The electron distribution param¬ 
eters are specified at tp = 0.5 and give an angular spread of Oq^ = = 2.0 and a transverse 

spread of = (Ty = 0.1. The electron trajectories are plotted in both x and y as T = 0 —)■ 1. 

The program produces an output graphic containing the results of the simulation, an example 
of which is shown in Figure 4.4. The output includes a list of the physical parameters used 
as well as a plot of the electron trajectories x{x) and y(T). In this simulation, = 0.5 and 
= TT, so the betatron oscillations complete only a half cycle over the length of the 
undulator and are focused at T = 0.5. 

4.4 Electron Trajectory Simulation Program 

The quadbeta program that I developed combines the functionality of quad and betaxy. It 
allows the simulation of electron trajectories through up to three quadrupole magnets as well as 
the undulator of an EEL. 

4.4.1 Program Structure 

The program begins by reading in the required parameters from either the command line or 
an input file. The program then initializes by seeding the random number generator, preparing 
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the graphical package to receive output data, and initializing variables. The input parameter 
window is also drawn at this time. 

Next, the electron trajectories are plotted in x. Sample electrons are generated with a Gaussian 
distribution in position and transverse velocity. The electron beam parameters are specified in 
the undulator, so the trajectories are first drawn there using a form of the analytical solution 
we determined in Equation (4.13). The electron trajectories are then plotted back through the 
quadrupoles, which are each modeled as a lens as we described in Section 4.1.2. The plot 
of x(t) is completed by drawing the quadrupoles as vertical purple lines and the start of the 
undulator as a vertical green line. The plot of y(T) is drawn in a similar fashion. 


Input Parameters 

As we mentioned above, the program can either read the required input parameters from the 
command line or an input file. Table 4.1 contains a list of the required input parameters with 
descriptions. 


Parameter 

(Opy 

"ffl, 'ff2, 'ffS 


Description 

Betatron oscillation frequency in undulator (x and y) 

Position in undulator at which beam parameters are specified 
Standard deviation in electron displacement from axis (x and y) 
Standard deviation in electron transverse velocity (x and y) 
Focal length of quadrupoles 1, 2, and 3 
Position of quadrupoles 1, 2, and 3 


Table 4.1: The input parameters required by quadbeta. These are either entered at the com¬ 
mand line or provided in an input file. All input parameters are dimensionless. 


The undulator is specified by the parameters and which give its betatron frequency in x 
and y. The electron beam parameters are specified at the position Tp with parameters describing 
the spread in positions (oy, Oy) and transverse velocities Oe^). 

The quadrupoles are described by their focal lengths (f i, f 2, f 3) and positions (Tfi, Tf2, Tfs). A 
quadrupole that is not used can be given a very large focal length (such as f 1 = 1 x 10^), which 
effectively removes any effect on the electron trajectories. 

Other parameters which can be controlled include the number of sample electrons to be simu¬ 
lated, the plot resolution of the output graphic, and the seed of the random number generator. 
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Simulation Output 

I have provided two sample outputs of the simulation: one modeling the effeet of a quadrupole 
doublet and another modeling the effeet of a quadrupole triplet. In both simulations I have 
speeified that the foeal point of the betatron oseillations oeeurs at t = 0.5, halfway through the 
undulator. This is often desired in praetiee—partieularly in an PEL oseillator, sinee the light is 
usually foeused near the eenter of the undulator by the optieal eavity mirrors. With the betatron 
oseillations also foeused there, the EEL interaetion is enhaneed, inereasing performanee. 

The quadrupole doublet in the first sample simulation is formed by two magnetie quadrupoles 
with equal strength arranged so that one foeuses in x and the other foeuses in y. A short drift 
spaee separates the quadrupoles, as was deseribed in Seetion 4.1.3. This results in a net foeusing 
effeet in both direetions, as ean be seen in Eigure 4.5. 

The simulation output in Eigure 4.6 shows how a quadrupole triplet is formed. The outer 
quadrupoles are foeusing in y and are half the strength of the inner quadrupole, whieh is fo¬ 
eusing in X. This eonfiguration also results in net foeusing in both direetions. 

Both quadrupole doublets and triplets are used in the beamlines of EEEs. Though a doublet ean 
provide net foeusing in both planes, a triplet provides equal foeusing in both planes, whieh is 
sometimes important [23]. 

4.5 Future Work 

This program has potential as a future design tool; it eould be used to assist in the plaeement of 
quadrupole magnets for optimization of the betatron oseillations in the undulator of an EEE, for 
example. It eould also be useful as a tutorial applieation for future students learning about EEEs. 

There are several ways that the program might be improved to inerease funetionality and po¬ 
tential applieations. Adding the simulation of dipoles would be a step toward extending the 
program to simulate a full beamline. Allowing more eomplex beam distributions would enable 
the study of effeets like beam halo (unwanted eleetrons well outside the primary eore of the 
beam, where they might strike beamline eomponents, potentially eausing damage and harmful 
radiation). The simulation eurrently models the quadrupoles as thin lenses and the fields in the 
undulator as having a hard edge. Simulating fields with finite spatial extent and soft edges would 
allow the study of phenomena like fringe field effeets for both the quadrupoles and the undu¬ 
lator. The program might also be extended to plot a beam envelope using the beam envelope 
equation in the quadrupoles and the undulator. 
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Electron Motion Through Quadrupoles and Undulator 
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x(T) 


-3 

3 


y (t:) 


-3 


Figure 4.5: Sample output of the eleetron trajeetory simulation program (quadbeta). A 
quadrupole doublet is formed with two quadrupoles separated by a small gap. Here f 2 = —f 3 
and the net foeusing effeet in both the x and y direetions is evident (the first quadrupole is ef- 
feetively “turned off” by letting f 1 = 1 x 10^). A list of input parameters is provided at the 
top, and below are plots of the eleetron trajeetories v(t) and The quadrupole positions are 
indieated with vertieal purple lines, and the start of the undulator is indieated with a green line. 


Input Parameters 

fi=le+09 £ 2 = 0.25 f3=-0.25 

Xfi=-1.6 Xf2=-1.55 Xf3=-1.5 

(Or =3.14 0 =0.1 On =2 

py y Qy 

=0.5 
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Electron Motion Through Quadrupoles and Undulator 
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Figure 4.6: Another sample output of the quadbeta program. A quadrupole triplet consists 
of a single quadrupole with a focal length / surrounded by two quadrupoles with focal lengths 
—2/. Just as in the quadrupole doublet case, properly selected focal lengths and quadrupole 
separation results in a net focusing effect in both directions. 


Input Parameters 

f,=0.4 f2=-0.2 £ 3 = 0.4 

Xfi=-1.6 Xf2=-1.55 Xf3=-1.5 

(Or =3.14 0 =0.1 On =2 
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CHAPTER 5: 
Conclusion 


The study I conducted on the effect of accelerator frequency on PEL performance shows that 
both amplifiers and oscillators perform well at lower accelerator frequencies. As the frequency 
is decreased, the average and peak electron beam currents are kept constant by increasing the 
length of the electron bunch, which is beneficial to both amplifiers and oscillators. Operating at 
a lower linac frequency also has the potential to reduce the overall cost, size, complexity, and 
power requirements of the EEL due to engineering advantages, including potential operation of 
an SRE linac at 4K. This study was conducted for infrared EEEs, but it could be extended to 
compare the applicability of these results at other wavelengths (e.g., terahertz, optical, or x-ray). 

Electron beam quality is critical to EEE performance, as the accelerator frequency study demon¬ 
strated. Beam quality can also affect how well the beam can be focused by quadrupole magnets 
throughout the EEE beamline. At the beginning of the undulator, quadrupoles can be used to 
focus the electron beam and ensure that the focal point of the betatron oscillations is near the 
center of the undulator to improve performance. The simulation I wrote to model electron trajec¬ 
tories in this process could be useful in the design of focusing systems for a new EEE—perhaps 
one that utilizes a lower linac frequency as suggested here. 

The EEE has the potential to provide the USN with the megawatt-class DE weapon it needs to 
deliver lethal energy at the speed of light in order to counter threats such as supersonic high- 
g-maneuvering ASMs. Significant advancements have been made since the USN first showed 
interest in DE weapons in the late 1960s, but many challenges remain. Continued research in 
key areas will be necessary before an EEE can be deployed on a ship. I hope that the work I 
have completed is a step towards making such a deployment a reality in the near future. 
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